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Abstract 

The preparation of neutron-optical phase gratings with hght-optical holography is reviewed. We compare the relevant 
concepts of i) Kogelnik's theory for Bragg diffraction of light by thick volume gratings, which can be used to analyze 
holographic gratings with both light and neutrons, and ii) the dynamical theory of neutron diffraction. Without going 
into mathematical detail, we intend to illuminate their correspondence. The findings are illustrated by analyzing data 
obtained from reconstruction of nanoparticle holographic gratings with both light and neutrons. 
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1. Introduction 



The photorefractive effect, i. e. the light-induced 
change of the light refractive-index in nonlinear opti- 
cal materials has been studied intensively since its dis- 
covery in 1966 lUj]. It was only in 1990 that it was 
realized by Rupp et al. [2] that its analog for neu- 
trons exists: The photo-neutronrefractive effect, i. e. the 
light-induced change of the neutron refractive-index of 
materials, has been investigated with cold neutrons in 
SANS experiments and used for neutron interferome- 
try |l3|-|6[] . Next to polymer-dispersed liquid crystals and 
deuterated polymethylmethacrylate a new material class 
seems promising for designing versatile neutron-optical 
devices such as beam-splitters and mirrors of high opti- 
cal quality: nanoparticle-polymer composites. In com- 
parison to holographic gratings recorded in the above- 
mentioned materials investigated earlier, the refractive- 
index modulation for neutrons of nanoparticle-polymer 
composites can be further adjusted by choosing the 
species of nanoparticles. Furthermore, these materials 
are resistant to volume shrinkage, which usually poses 
a problem in polymerization processes |l7[]. 
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Figure 1 : Recording of liolograpliic gratings with laser light, d and A 
are the mechanical thickness and spacing of the grating, respectively. 



2. Recording of holographic gratings with light 

In Fig.[T] the experimental setup for recording holo- 
graphic gratings in photosensitive materials is sketched. 
Two coherent plane s-polarized light waves of ampli- 
tude Ao/2 and wavelength Ai - In/ki in free space, 
whose y-component can be written as 

(As = Ao/2 exp [iki (z cos 6^"^ + x sin e*""*)] , 
iffR = Ao/2 exp [iki (z cos 0<''> - x sin e^^'^^j , 

are brought to interference within the material. Each 
beam encloses the (external) incidence angle 0^'"' with 
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the sample surface normal. The resulting light pat- 
tern exhibits a sinusoidal modulation I - \^s + "ArP = 
A2/2[1 +co5{Klx)], where Kl = 27r/A = 2fez.sin0<"'. 
For a given material, illumination time and intensity, the 
light-induced refractive-index change may be written as 
An/,(jc) = Arv^ + Arv^ cos,{Kix). In the present case, 
the photosensitive material consists of a nanoparticle- 
monomer-photoinitiator mixture IlZHsI] . The photoini- 
tiator triggers polymerization upon illumination with 
light. In regions with high intensity, monomers are con- 
sumed by the polymerization reaction. More monomers 
diffuse from dark to bright regions. Therefore, nanopar- 
ticles (e. g. Si02) migrate to dark regions. 



3. Diffraction from holographic gratings 

The recorded holographic gratings can be used for 
diffraction of light or neutrons. In 1969, Kogelnik 
OlOll formulated a two-wave coupling theory for Bragg 
diffraction of monochromatic light by thick sinusoidal 
holographic gratings starting from the Helmholtz equa- 
tion for the complex electric field amplitude inside the 
grating. The approach is valid for reflection and trans- 
mission holograms and also at large diffraction efficien- 
cies rjL = IdI{Id + If), with the intensity of the Bragg- 
diffracted beam Iq and the intensity of the beam in for- 
ward direction If. It includes phase as well as absorp- 
tion gratings. Moreover, also slanted geometries - in 
contrast to the symmetric situation in Fig.lT] in which 
no refractive-index modulation occurs in z-direction - 
can be treated easily. In particular, the theory assumes 

1 . incidence near the Bragg angle, 

2. only two coupled waves propagating inside the ma- 
terial, 

3. slow energy variation in space of the two coupled 
waves compared to their wavelength (be it due to 
energy exchange or absorption), 

4. that second derivatives of wave amplitudes can be 
neglected (this is a consequence of the previous as- 
sumption), 

5. that outgoing waves are equal to the waves inside 
the material. In other words, boundaries between 
vacuum or air (light refractive-index equal or close 
to 1 ) and the material {mean light refractive-index 
larger than 1) are neglected. The mean refractive- 
index is assumed to be constant in space. Only the 
modulation of the refractive-index A«^ is taken 
into account. This is again a consequence of the 
previous assumption, as is discussed in llllll . 



These assumptions lead to two first-order differential 
equations (coupled wave equations), either containing 
both wave amplitudes. 

3.1. Light 

Solving the coupled wave equations mentioned above 
for the special case of lossless dielectric holographic 
unslanted gratings lIlOll . the light diffraction efficiency 
reads as 



riL = v^sinc^ 



(V^l^)^ 



(1) 



in transmission geometry (Laue case), with 



Vl = K 



COS 6 Ai cos 6 



nAn^l^d ndids-O) 



A 



Here, A«^'\ d^B = arcsin [/tz,/(2Anz.)] and 6 = 

arcsin(sin0'^''V«L) are the amplitude of the light- 
induced sinusoidal refractive-index modulation for 
light, the Bragg angle and the actual angle of incidence 
(both as measured inside the material), respectively, ni 
is the mean refractive-index of the material for light 
including A«^ . k is the coupling constant giving the 
strength of coupling between the two waves that propa- 
gate through the material. A«^ arises from the modu- 
lation of the dielectric permittivity e due to inhomoge- 
neous distribution of nanoparticles and polymer. 

3.2. Neutrons 

From the time-ir idep endent one-particle Schrodinger 
equation (see e. g. 11211 ') 



_^V2 + V{f) 
2m 



^(^ = E,^(i), 



one obtains a Helmholtz-type wave equation 
[V^ + K^if)] i//(i^ = 



IIJ, with K{r) = ^l2m[E- V{r)\lh, the neutron wave 
number inside the material. E - ffk^/{2m) is the en- 
ergy of the free neutron. As usual, one defines the neu- 
tron refractive-index of the material as n - K/k, with the 
wave number k - 2n/A of the incident neutron beam. 
The neutron-optical potential in nanoparticle-polymer 
holographic gratings is 



V = 27Th^bc [p + Ap(x)] /m, 



(2) 



with the neutron mass m, the coherent scattering length 
be for a particular nanoparticle-polymer mixture and 
the number density p. The photo-neutronrefractive ef- 
fect in nanoparticle-polymer composites arises from the 



modulation of the neutron-optical potential due to in- 
homogeneous distribution of nanoparticles resulting in 
Ap(x) + 0. The neutron refractive-index becomes 






1 



??hc \p + Ap(x)] 
In 



and therefore the neutron diffraction efficiency r\ can be 
written in a similar form as in Eq. ([TJ, with vl and ^z, 
replaced by 

Adbc^p ^ nd{0B-B) 



^ 



(3) 



2 cos A 

respectively. Here, the Bragg angle is Ob - 
arcsin[/l/(2A)]. From the above equations, it is ex- 
pected that the following parameters can be adjusted for 
fabrication of neutron-optical elements: First, the grat- 
ing thickness d can be freely chosen up to a material- 
dependent value at which extinction inhibits penetra- 
tion of the material with light and therefore hologram 
recording. Second, the grating spacing A is determined 
by 0^''' for given A^. Unfortunately, A„„„ = Al/2 - 
which is at least in the range of lOOnm for laser wave- 
lengths - is the lower limit for given A^. This results 
in very small Bragg angles for neutrons. A possible 
solution to overcome this limit is triggering nonlinear 
processes in the material that generate higher harmonics 
of the refractive-index modulation, as has already been 
demonstrated in 1 14). Third, Ap can be adjusted by the 
illumination time and/or intensity. Finally, the coher- 
ent scattering length b^ is determined by the particular 
choice of nanoparticles and their isotopes. 

4. Correspondence to the dynamical theory of neu- 
tron diffraction 

Dynamical diffraction theory for scattering of neu- 
trons (see e. g. lll2l Il3l Il5ll ) in crystals is based upon 
solving the one-body Schrodinger equation outside and 
inside the material, the latter being represented by a pe- 
riodic optical potential. It is assumed that the Born ap- 
proximation is invalid, i. e. multiple scattering effects 
are taken into account. Furthermore, only two waves 
- forward diffracted and Bragg-reflected - are assumed 
to have considerable amplitudes. Employing a superpo- 
sition of Bloch modes as ansatz and taking the bound- 
ary conditions into account, one finds two values for the 
refractive-index for each of the two internal waves and 
so the coherent wave is a superposition of four modes. 
One can calculate the reflectivity R of a thick crystal in 
the symmetric Laue case as: 



R(x,y) 



'■{y^/^TT} 
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Figure 2: ± 1 ^' order rocking curves (left/riglit) of a nanoparticle liolo- 
grapliic grating witli 20 voI% SiOo nanoparticles, d ^ SQ/ivn and 
A = 1 /im measured witli laser light of wavelength Al = 458 nm. 



with the dimensionless parameters 

7rVsin(20B) nd Ad\F\ 

(x and y should not be confused with the spatial co- 
ordinates X and y in Fig.[Tli. Here, V and \F\ are 
the volume of the primitive unit-cell and the unit-cell 
structure factor for a particular reflection plane, respec- 
tively. The quantity A is called extinction length or Pen- 
dellosung period. It determines the oscillation period of 
the internal amplitudes of forward diffracted and Bragg- 
reflected waves along the direction parallel to the sur- 
face normal (z-direction in Fig.[T]). Pendellosung inter- 
ference fringes for Laue reflections from silicon crys- 
tals have been demonstrated by variation of the incident 
wavelength Ill6tl . The effect is also present in oscilla- 
tions of the integrated reflectivity, that have been shown 
iniH. 

Comparing Eqs. ([T]) and (|4| one can see the following 
correspondence between the parameters v, ^ used in the 
Kogelnik theory and y, x used in the dynamical diffrac- 
tion theory for neutrons: 

y ^ V 
X ^ i 

V 

Calculating an analogue A^ to A via the relation v = 
nd/AK, one finds 

2jtcos0b 



AbcAp 



(5) 



5. An example 



(4) 



Two glass microscope slides are glued together en- 
closing a fixed gap provided by a roughly 50 /im thick 
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Figure 3: Diffraction pattern of the same sample as analyzed in Fig.|2] 
measured by analyzer crystal rotation at a neutron wavelength of about 
1.9 A. 
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Figure 4: Rocking curve of the same sample with analyzer crystal at 
-1^' order peak position (c. f. Fig.[3j. 



spacer-foil. The liquid monomer-nanoparticle mixture 
is sucked into the gap by capillary force. After record- 
ing a hologram in the material as described in Sec.|2] a 
laser beam of wavelength Al = 458 nm is diffracted by 
the holographic grating with 20 vol% of SiOa nanopar- 
ticles of roughly lOnm diameter, grating thickness d - 
50 /im and grating spacing A - I jjm. By rotating the 
sample through the Bragg-positions for -P' and -i-P' 
diffraction orders, the rocking curves shown in Fig.|2]are 
observed. A fit according to Eq. ^ (solid lines) yields 
d = 50.60(5) /im, while the first-order refractive-index 
modulation for light is found to be An['^ = 5.007(9) X 
10"^. Absorption is negligible at this wavelength. 

The sample has also been analyzed with thermal neu- 
trons using the USANS-option of the interferometer 
instrument S 18 of the Institute Laue Langevin (ILL), 
Grenoble [18]. This double perfect crystal diffractome- 
ter setup accesses scattering angles down to the //rad 
range, corresponding to /im structures in real space. The 
wavelength of the incident beam was A ^ L9 A, roughly 
four orders of magnitudes less than the grating spacing 
of the sample. By rotation of the analyzer crystal at nor- 
mal incidence of the beam to the sample, the diffraction 
pattern as shown in Fig. [3] is measured. Fixing the an- 
alyzer crystal position at the -P' order peak (left peak 
in Fig.O and rotating the sample through the -P' or- 
der Bragg-angle, the rocking curve as shown in Fig.|4] 
is observed. The solid line is a fit according to Eq. ([TJ 
including Eqs. (O, which describes well the measured 
data. The fit yields the parameter values d - 49.1(6) //m 
and bcAp - 4.36(5) /zm"^. The Bragg angle at this 



wavelength is Og - 6 x 10"^ °. One can calculate the 
Pendellosung period using Eq. (|5]l to find Afc - 7 mm, 
a value that is about two to three orders of magnitude 
larger than for the (220) reflection of a silicon crystal 
slab. The Pendellosung oscillations are described by 



TliOB)-- 



sin 



^'^^^■d]»Ries)^sini^^.d 
2 cos 6b j \ V cos 6b 



within the frame of the two corresponding theories com- 
pared in this paper From the above expressions it im- 
mediately becomes clear that this enormous value is 
mainly due to bcAp being very small in holographic 
nanoparticle-polymer gratings. 

6. Discussion and concluding remarks 

It is important to note that three diffraction orders 
are clearly visible in Fig. [3] so obviously three waves 
instead of two - as was assumed in Sec. [3]- are cou- 
pled propagating through the sample. Therefore, the 
grating can be considered to be a 'thin' grating at this 
wavelength, i. e. diffraction occurs in the multi-wave 
coupling- or Raman-Nath-regime and not in the two- 
wave coupling- or Bragg-regime (see e. g. [191 ). How- 
ever, since the samples' diffraction efficiency for ther- 
mal neutrons is very small, the intensity that contributes 
to the additional (third) peak is negligible and two-wave 
coupling theory can be applied to analyze the measured 
rocking curve in Fig.|4] For larger values of 77, rigorous 
coupled wave analysis |20] must be employed if more 
than one reflected beam is observed. 



It is well known that for thermal neutron diffraction 
by a thick silicon crystal in the Laue case, the reflec- 
tivity R in Eq.|4] oscillates rapidly as a function of the 
parameter x, that reflects the deviation from the Bragg- 
condition 9 -9b for a given crystal and wavelength. This 
is due to the high values of y, that are typically more 
than 10;7r. As a consequence of limited resolution, the 
measured reflectivity can be at most R - 0.5 at Bragg- 
angle. One of the major advantages of holographic grat- 
ings over perfect crystals is that such limitations can be 
avoided. The parameter v - corresponding to j/ - that 
determines the frequency of oscillation of R depending 
on 0- 9b can be tuned by careful choice of bed and Ap 
and is typically smaller than n. As a consequence, the 
diffraction efficiency rj or reflectivity R can well exceed 
the Umit of 0.5 and can, in principle, be equal to 1 at 
V - 71 12. For instance, one can implement mirrors in- 
stead of the usual beam-splitters for the second grating 
of triple-Laue neutron interferometers. Used that way, 
mirrors save intensity that is lost in perfect-crystal neu- 
tron interferometry. Mirrors also prevent contribution to 
the background by parasitic beams whose signal cannot 
be separated from the signal of the interfering beams in 
the case of small Bragg-angles |5t]. 
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